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We present a master formula, with applicability conditions, which allows 
us to automate the resummation of infrared and collinear logarithms 
appearing in distributions of jet observables in QCD at next-to-leading 
logarithmic accuracy. 

1 Jet observables: fixed order and resummation 

Jet observables, event shapes and jet rates, revealed themselves as one of the 
richest laboratories to explore QCD. Being infrared and collinear safe (IRC), 
they can be predicted with perturbative (PT) techniques, but their high sensi- 
tivity to low energy emissions allows us to investigate the fairly unknown non- 
perturbative regime. Most discriminatory studies make use of distributions. 
In integrated distributions one requires that the value of the observable 
V({p}, ki, ■ ■ ■ , k n ) - a function of all secondary final state momenta fej and of 
the Born momenta after recoil from the emissions {p} - be less than a fixed 
value v , , 



The inclusive phase space region, where v — O (1), is dominated by events with 
hard jets, which move the value of the observable far away from its Born value, 
v = 0. These events can be described with fixed order PT expansions, however 
they are quite rare, every additional jet being suppressed by an additional 
factor of cv s . More common events are characterised by a large number of soft- 
collinear emissions which modify only slightly the value of the observable from 
its Born value, so that » < 1. Here fixed order predictions fail since every 
power of a s is accompanied by up to two large logarithms L = ln(l/w) of the 
value of the observable. A reorganization of the PT expansion is then needed 
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in order to resum all leading (LL, exp{a"L" +1 }) and ncxt-to-leading (NLL, 
exp{a™£™}) logarithmic terms. 

In the past few years the analytical resummation for a variety of observ- 
ables have been presented in e + e~-collision[l, 2] and DIS[3]. The matching 
of resummed predictions with fixed order results allowed tests of QCD, mea- 
sures of the coupling constant and studies of non-perturbative corrections [4]. 
However, the need for a separate analytical calculation for every observable 
has limited the experimental use of resummed predictions. When dealing with 
multi-jet observables, analytical calculations become quite unfeasible, involv- 
ing many integral transforms in order to write the distribution in a factorized 
form [5]. Also in some cases it turns out not to be possible to resum the ob- 
servable analytically [6]. 

We present then here a general approach to resummation based on a prelim- 
inary automated analysis of the observable, to establish its relevant properties 
with respect to soft-collinear emissions; in a subsequent step this information 
is used as an input of a general master formula. 

2 Automated resummation 

2.1 Analysis and applicability conditions 

We start by considering a Born event consisting of n hard partons ('legs'), n, of 
which are incoming, with momenta p\ . . .p n . We resum (n+l)-jet observables 
in the n-jet limit. The (positive defined) observable should then 

1. vanish smoothly as a single extra (n+1) parton of momentum k is made 
soft and collinear to a leg £, with the functional dependence 



Here Q is a hard scale of the process and the secondary emission k is de- 
fined in terms of its transverse momentum k t and rapidity r] with respect 
to leg £, and where relevant, by an azimuthal angle <j> relative to a Born 
event plane. By requiring the functional form (2) (in practise, almost 
always valid), the problem of analyzing the observable reduces in part 
to identifying, for each leg £, the coefficients ae, be, de and the function 
ge{4>). IRC safety demands ae > max{0, — be}. 
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2. be recursively IRC safe, i. e. given an ensemble of arbitrarily soft and 
collinear emissions, the addition of a relatively much softer or more 
collincar emission should not significantly alter the value of the observ- 
able, condition required for exponentiation of the leading logarithms. Ob- 
servables like jet-rates in the Jade algorithm are then excluded, and many 
other examples of non-exponentiating IRC safe obscrvables exist [7]. 

3. be continuously global[3] — this means that for a single soft emission, 
the observable's parametric dependence on the emission's transverse mo- 
mentum (with respect to the nearest leg) should be independent of the 
emission direction. In practice this is perhaps the most restrictive of the 
conditions. It implies ai = a 2 = . . . = a n = a. 

To enable our computer program CAESAR (Computer Automated Expert 
Semi- Analytical Rcsummation) to establish these properties with the desired 
degree of reliability and precision, we have found it useful to make use of 
multiple-precision arithmetic [8] . 

2.2 Master formula 

Given the above conditions, one can derive the following NLL master rcsum- 
mation formula for the distribution E(u)[7] 

n 

InS(«) = -£C/ 
i=i 

+ E ln tf) +lnS {T{L l a ^ + ^(CiK , . . . , C n r' n ), (3) 

where Ci is the color factor associated with Born leg £, and Et is its energy, 
B( is —3/4 for quarks and — (11CU — 4Tr?i/)/ '(12CU) for gluons, Inde = \nd? + 
J *^hj/(^), and for incoming legs, fi are the parton densities. Wc note 
that (3) is independent of the frame in which one determines the di and (to 
NLL accuracy) of the choice of hard scale Q. 

The functions rt(L), which contain all the LL (and some NLL) terms, are 
rv — Q 2 dk 2 as(kt) / h \ tj rQ 2 d kfa s (k t ) Q 



r e (L)+r' e (L) find/ -6/ In 



2Et 

Q 



+ B e T 



L 



be 
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here a s , in the Bremsstrahlung scheme [2], runs at two-loop order. r' e = dire 
and T(L) are relevant only at NLL level 

T(L)=[ Q2 d 4^±. (5) 

J e -2LQ2 k t TT 

The process dependence associated with large-angle soft radiation is de- 
scribed by S(T(L/a)), whose form depends on the number n of legs. For n = 3 



\nS(t) = -t 



qgQq'g or< , Qqq' 



W <?<?' 



(6) 



where Q 2 ab = 1p a .pb and q, q' and g denote the (anti)-quarks and gluon. The 
simpler case n = 2 can be read from (6) setting Ca — 0. The n = 2, 3 formulae 
apply then to e+e~, DIS and Drell-Yan production, while a process such as 
gg — > Higgs + (/ involves simply different color factors. The (more involved) 
n = 4 case needed to describe hadronic dijet production can be found in [7]. 

We examine now the factor T . Without it, eq. (3) corresponds essentially 
to the probability of vetoing all emissions k with V({p},k) > v. However 
at NLL subtle effects enter: a simple veto on single emissions turns out to be 
insufficient, since events might have k±, . . . , k m ) > v though all emissions 

separately had fcj) < v, or vice versa. This multiple emission effect, 

encoded in T, is connected to how all secondary emission coherently determine 
the value of the observable and can be computed in a general way in [6]. 

While this talk concentrates on the method, new results obtained with it 
were presented at this conference in[9] (and other results can be found in[7]). 



3 Final remarks 



This project started with the development of an algorithm to numerically com- 
pute the non-trivial NLL terms associated with multiple emissions. It was ini- 
tially applied by hand to three new observables [6] . Progress since then includes 
the derivation of a general master formula, together with a precise, automati- 
cally verifiable list of conditions on the observable for the resummation to be 
valid at NLL and the full numerical implementation of this. Though no human 
intervention is needed, results have the quality of analytic NLL predictions, so 
that any hadronisation model can be applied, studies of renormalization and 
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factorization scale dependence can be carried out and a matching with fixed 
order is feasible. The most import results obtained up to now are the first re- 
summations in hadronic dijet production, for a variety of observable at a time. 
We now aim at automating the matching with fixed order results[10], this will 
open up the possibility to carry out a vast amount of phenomenological studies. 
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